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Abstract 



Velocity measurements of wind blowing near the North Sea border of North- 
ern Germany and velocity measurements under local isotropic conditions of a 
turbulent wake behind a cylinder are compared. It is shown that wind gusts - 
measured by means of velocity increments - do show similar statistics to the 
laboratory data, if they are conditioned on an averaged wind speed value. 
Clear differences between the laboratory data and the atmospheric wind veloc- 
ity measurement are found for the waiting time statistics between successive 
gusts above a certain threshold of interest. 



1 Introduction 

The occurances of gusts represent a great problem affecting a wide field of technical 
applications as for instance the construction of Wind Energy Converters |I| as well 
as the construction of high buildings in general [fj] . But also regarding the landing- 
and take-off-operations at airports the importance of the statistics of gusts becomes 
evident ||. Nevertheless gusts are not uniquely denned. In j| or || gusts are taken 
as an additional structure besides the permanent, fluctuating turbulence u(t) around 
the mean wind u(t) of the wind velocity field. More often wind gusts are charac- 
terized by extrem velocity differences during short time intervals ranging typically 
between some seconds up to a minute. But these characterisations are still vague. 
In order to examine wind gusts in a statistical way we use a data set of a turbulent 
wind velocity field recorded near the coastline of the Northsea at the FH Ostfries- 
land in Emden (Germany). The velocity was measured by means of an ultrasonic 
anemometer at 20 m height. The sampling frequency was 4 Hz. The measuring 
period took about one year (1997-1998). After careful investigation of the quality of 
the data we examine a representative 275-hour-excerpt of October 1997 and focus 
on the velocity component in direction of the mean wind (See also for a first brief 
report on the analysis of these data ||). 

Our main interest is to explore how atmospheric turbulence is related to the well 
known local isotropic and stationary laboratory one. Therefore we use a simple 
measure for gusts and compare both velocity fields by means of their probabilistic 
features. The laboratory data we use here was recorded in a wind tunnel behind 
a cylinder with a Reynoldsnumber of Re ~ 30000 (see [0]). In a second step we 
examine the waiting time distributions of successive wind gusts to resolve also their 
time structure. 
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2 Probabilistic description of wind gusts 

The wind field is known to exhibit a high degree of turbulence with Reynoldsnumbers 
Re of about 10 7 . The wind velocity U(t) is commonly (see Q) expressed as the sum 
of the mean velocity u(t) and fluctuations u(t) around it. For the mean velocity u(t) 
a ten minute average is used: 

U = u + u . (1) 

Here we concentrate on the velocity components in the mean velocity direction. 
The greater the fluctuation values u the more turbulent the wind field becomes. 
In windenergy research this is often expressed by the turbulence degree ti which is 
defined as the standard deviation a in relation to the mean velocity u: 

U = \. (2) 

Nevertheless the value of ti does not contain any dynamical or time-resolved infor- 
mation about the fluctuation field itself. 



v(t) A 
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Figure 1: Idealized representation of a wind gust. 



To achieve a deeper understanding of wind gusts - as a result of the fluctuating wind 
field - we investigate in how far wind gusts are related to the well known features 
of small scale turbulence. To this end we will perform an analogous analysis for 
our wind data and data from a turbulent wake flow behind a cylinder in a wind 
tunnel. As a natural and simple measure of wind gusts we use the statistics of 
velocity increments 5u T of the fluctuations: 

5u T := u{t + r) - u{t) . (3) 

The increments directly measure the velocity difference after a characteristic time r 
(illustrated in Fig.|). So a high increment exceeding a certain threshold S (Su T > S) 
can be defined as a gust. 

For a statistical analysis we are interested in how frequent a certain increment value 
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occurs and whether this frequency depends on r. Therefore we calculate the prob- 
ability density functions (pdfs) P(5u T ) of the increments. In Fig. |2] the pdfs for 5 
different values of r are shown. These distributions are all characterized by marked 
fat tails and a peak around the mean value. Such pdfs are called intermittent and 
differ extremely from a Gaussian distribution that is commonly considered to be the 
suitable distribution for continuous random processes. 



P(5u x ) / a.u. 
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Figure 2: The pdfs of the unconditioned windincrements for r being 0.008 T, 0.03 T , 
0.2 T, 0.95 T and 61 T (full symbols from the top to the bottom) are drawn in. 
The pdfs are shifted against each other for a clearer presentation and the respective 
fitfunctions according to eq. (fjj are shown as solid curves. 



A normal distribution is uniquely defined by its mean value /i and its standard 
deviation a. Thus every distribution can be compared to a normal distribution in a 
quantitative way. In Fig. |3| we compare one of the measured pdfs (r = 4 s) with a 
Gaussian distribution with the same a. In this presentation the different behaviour 
of the tails of both distributions becomes evident. Note that the high increments 
of the wind pdfs - located in the tails - correspond to strong gusts. For instance 
the value of 5u T = 7 • cr corresponds to a velocity ascending of 5.6 rajs during 4 s). 
As shown in Fig. ^| the measured probability of our wind data is about 10 6 times 
higher than for a Gaussian distribution with the same standard deviation! This case 
is represented by the arrow in Fig. [| The value 10 6 - for instance - means that a 
certain gust which is observed about five times a day should be observed just once 
in 500 years if the distribution were a Gaussian instead of the observed intermittent 
one. 

But intermittent distributions seem to appear quite often in natural or economical 
systems like in earthquake- (0), foreign exchange market- (||) or even in some 
traffic-statistics ([|T0|j). 
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a-P(5u x /a) 




5u T / a 



Figure 3: The distribution of the windincrements for r = 4s is represented by the 
squares, a Gaussian distribution with the same standard deviation a by the solid line 
(parabola due to the semilogarithmic presentation) . Both distributions are normal- 
ized with a = 0.8 rajs. 



What kind of statistics do we get in the case of local isotropic and stationary labora- 
tory experiments? The typical probability density functions in laboratory turbulence 
- as shown in Fig. [| - change from intermittent ones for small values of r to rather 
Gaussian shaped distributions with increasing r. When r rs T with T being the 
correlation time (integral time) the approach to a normal distribution is marked: 



T = J R(r)dT . (4) 
o 

R(t) is the correlation function of the fluctuations. In our case the correlation time 



of the atmospheric flow is 34 s and 6 ms for the laboratory data [D]] . 

Note that for the pdfs of the atmospheric velocity field this characteristic change of 

shape, even for r-values higher than T (as shown in Fig. 0) is not observed. 
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Figure 4: In a) the symbols represent the pdfs P(5u T ) of the laboratory increments 
for different values of r. From the top to the bottom r takes the values: 0.005 T, 
0.02 T , 0.17 T, 0.67 T and 1.35 T. In frj the conditioned pdfs are presented, here r 
is 0.008 T, 0.03 T , 0.2 T, 0.95 T and 1.9 T. T/ie mean wind interval on which the 
increments are conditioned is [4.5 ; 5.6] m/s. In both cases the solid lines are the 
corresponding fit-functions according to eq. (|3p. The distributions and their fits are 
shifted against each other for a clearer presentation. 



As already mentioned a fundamental difference between atmospheric and laboratory 
turbulence is that the latter is stationary. In laboratory experiments one usually 
deals with a fixed value and direction of the mean wind speed u, which is obviously 
never the case for atmospheric wind fields. Therefore in a second step we calculate 
the pdfs of the atmospheric increments only for certain mean velocity intervals. That 
means that only those increments are taken into account with u ranging in a narrow 
velocity interval with a width of typically 1 m/s. These conditioned pdfs P(5u T \u) 
show a similar qualitative change of shape like those of the laboratory experiment^ 
which is illustrated in Fig. |] b). 

To quantify this similarity we use a well established fit by an empirical explicite 
function for the pdf. This formula was derived in [E| on the basis of Kolmogorov's 



1 Only for very small values of u (u < 1 m/s) this change of shape is not observed. 
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understanding of a turbulent cascade: 

JL °f 5u 2 ln 2 {s/s ) d{lns) 
P{5u T ) = J exp{-^) • exp( ^~)— ■ (5) 

o T 

In Fig. |2| and ^ these fitfunctions are represented by the solid lines. A 2 , is the 
fundamental parameter (formparameter) in equation ([5]) and determines the shape 
of the probability distribution. As it can easily be seen equation fl5|) reduces to a 
Gaussian distribution if X 2 is zero: 



lim P(5u T ) = j==exp{— t^-j) • (6) 



1 , Su 2 

On the other hand the more X 2 increases the more intermittent the distributions 
become. In this way the parameter X 2 may serve to compare the pdfs with each 
other in a more quantitative way. In Fig. |^ the evolution of the formparameter as 
a function of the increment distance r is shown. 

Other laboratory measurements (see [12, pj) of X 2 have shown evidence that it 



saturates approximately at 0.2. As shown in Fig. the formparameter of the con- 
ditioned wind increments as well as of the laboratory ones is approximately 0.2 for 
small r— values. Furthermore it tends to zero with increasing r. None of these two 
features is observed in the case of the unconditioned increments, X 2 is rather inde- 
pendent from t with a value of about 0.7. A constant behaviour of A^ means that 
the form of the pdfs remain unchanged, while its variance may change. Together 
with the power law scaling of < 5u 2 ~ r 2 / 3 > we get the unusual result, that for pdfs 
with constant intermittent shapes the scaling of the n-th moments < 5m™ ~ r^ n > 
follow the proposed law of Kolmogorov 41|TJ]], £ n = n/3, i.e. without intermittency 



corrections for the exponents. 

Thus we have shown that the anomalous statistics of wind fluctuations on discrete 
time intervals - which are obviously related to wind gusts - can be reduced to the 
well known intermittent (anomalous) statistics of local isotropic turbulence. This 



result deviates from results of wind data reported in [15], where it is claimed that 
their unconditioned wind pdfs behave like stationary ones from laboratory measure- 
ments. Thus on the basis of our stochastic analysis we claim that wind gusts are no 
independent structures but are a part of small scale turbulence. 
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Figure 5: Aj as a function of the increment distance r is shown for (from top to 
bottom) the unconditioned, the conditioned and the laboratory increments. For the 
last r = 1 means 10~ 5 s for the two former 0.25 s For a clearer visualisation a 
double-logarithmic presentation is used. 



3 Waiting time distribution 

So far we have shown how the atmospheric turbulence is related to the laboratory 
one in a statistical way. This probabalistic ansatz describes the frequency with 
which certain gusts occur but it is not clear how they are distributed in time. In 
this sense we now examine the inter-event-times between successive wind gusts. 
The marked fat tail behaviour of the unconditioned pdfs - as illustrated in Fig. |2| 
a) - points at an interesting effect. In |16| the equivalence between the divergence 
of the moments < x q > and the hyperbolic (intermittent) form of pdfs which leads 
to a power law behaviour of the probability distribution is emphasized: 

p{x>S)(xS- q ,S»1 . (7) 

A famous example of such a natural power law behaviour is the Gutenberg- Richter- 
law |H| that describes the frequency N of earthquakes with a magnitude being 
greater than a certain threshold M (magnitude): 

log(N(E)) = a-bM 
<&N(E) oc E~ h . (8) 

N(E) describes the number of occurances of those earthquakes with an energy out- 
put higher than a threshold Eq. The magnitude M is proportional to the logarithm 
of E. The same result is found for the energy flux distribution of a wind field ||. 
But also the waiting time distribution of fore- and after shocks obey a power-law, 
what is known as the Omori-law [|T8"| . 
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In this sense we now examine the waiting time distribution of wind gusts. Therefore 
we refer to the idealized gust represented in Fig. [TJ choosing different thresholds A 
and different increment distances r (see eq. [3]). Always when the condition 5u T > A 
is fulfilled a gust event is registered. To avoid that one event is counted several 
times we use the condition that the distance between two successive events AT is 
at least r. Due to the cumulation of occurances for small time distances we choose 
a logarithmic time. 

The distributions for A = 4.0 m/s and r = 0.3 T (about ten seconds) and for 
A = 1.5 m/s and r = 2.0 T are shown in Fig. |6| a) respectively || b). It can easily 
be seen that the exponential fit in Fig. |6| corresponds to a power law concerning the 
waiting times AT. Only the exponent seems to depend on A and r. 
Considering 

x = ln(AT) (9) 

and the normalisation 



1 = / P(x)dx = J P(x(At))^-d(AT) , (10) 

^ ' 

the used fitfunction 



P(x) oc e~ cx (11) 
can also be written as a function of AT which leads to a power-law: 

P(AT) oc k ■ e^ AT ^~ c ■ AT' 1 

oc k-(AT)- - 1 . (12) 



This power law behaviour of the waiting time distributions is only observed for the 
atmospheric wind data and not for the stationary laboratory one. 
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Figure 6: The filled symbols illustrate the distribution of the logarithm x of the 
waiting times AT between successive gusts. Additionally the fitfunctions are drawn 
in (solid lines), x = 8 for instance corresponds to a time distance of AT m 1000 s. 
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4 Summary 

On the basis of well denned velocity increments an analogous analysis of measured 
wind data and measured data from a turbulent wake was performed. The statistics 
of velocity increments, as a related statistics to the occurance frequency of wind 
gusts, showed that they are highly intermittent. These anomalous (not Gaussian 
distributed) statistics explain an increased high probability of finding strong gusts. 
This could be set in analogy with turbulence measurements of an idealized, local 
isotropic laboratory flow if a proper condition on a mean wind speed was performed. 
So far the statistics of wind gusts can be interpreted as a superposition of idealized 
packets of turbulence at different Reynoldsnumbers. 

As a further statistical feature of wind gusts we have investigated the waiting times 
between successive gusts exceeding a certain strength. Here we find power-law- 
statistics (fractal statistics) - similar to earthquake statistics - that can not be 
reproduced in laboratory measurements. 

To conclude we have shown two important aspects of wind gusts. The overall oc- 
curance statistics could be set into analogy to the anomalous statistics of velocity 
increments in local isotropic turbulence. The time structure of successive gust events 
displays fractal behaviour. 
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